DESY-97-060 



Measurement of Elastic J /ifj Photoproduction at 

HERA 



ZEUS Collaboration 



Abstract 

The reaction 7 p — > J /ip p has been studied in ep interactions using the ZEUS 
detector at HERA. The cross section for elastic J/il) photoproduction has been 
measured as a function of the photon-proton centre of mass energy W in the range 
40 < 1^ < 140 GeV at a median photon virtuality of 5 x 10"^ GeV^. The 
photoproduction cross section, o'^p_>jy^p, is observed to rise steeply with W . A 
fit to the data presented in this paper to determine the parameter 5 in the form 
o"^p^j/^,p oc yields the value 5 = 0.92ib0.14ib0.10. The differential cross section 
da /d\t\ is presented over the range \t\ < 1.0 GeV^ where t is the square of the four- 
momentum exchanged at the proton vertex, da /d\t\ falls exponentially with a slope 
parameter of 4.6 it 0.4lQ'g GeV~'^. The measured decay angular distributions are 
consistent with s-channel helicity conservation. 
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1 Introduction 



This paper reports new data on the photoproduction of the J/ ip meson using the ZEUS 
detector at HERA. It is part of our continuing study of vector meson {p,u,(f), J/ip) pro- 
duction in both the photoproduction |l|, |[ and the deep inelastic scattering regimes 
^. Previous results have established a weak dependence on the photon-proton centre 
of mass energy, W, of the vector meson photoproduction cross sections (a oc with 
5 ~ 0.2 — 0.3) if there is no hard scale in the process, as expected from soft diffraction. 
By contrast, the cross sections for elastic p and (p production in deep inelastic scatter- 
ing at 5<Q^<20 GeV^ exhibit a stronger W dependence (a oc with 6 ~ 0.3 — 0.6) 
where sets the hard scale. In the photoproduction of the J/ip meson the mass of the 
J/ip itself provides the hard scale and the cross section exhibits a strong W dependence 
(cr oc with 5 ~ 1). The total virtual photon-proton cross section 0, |] also exhibits 
a change in energy dependence as increases beyond ~ 1 GeV^. Overall, the data 
illuminate the transition from the soft, non-perturbative regime to the kinematic region 
where perturbative descriptions become applicable. 

J/tp photoproduction has been measured as a function of W from threshold to W ^ 
20 GeV in fixed target experiments g and extended to W ^ 140 GeV at HERA 

|T2|, . A review of the low energy experimental results can be found in reference [|1^ 



In this paper we extend our earlier study of elastic J/ip photoproduction to include 
the determination of the differential cross section da/d\t\ and the angular distributions of 
the decay leptons. In addition, the six-fold increase in the size of the data sample allows 
us to determine the parameter 6 from the data presented here alone. 

The J/ ip was detected via its leptonic (electron pair and muon pair) decay modes in the 
kinematic range < W < 140 GeV. After a brief description of the ZEUS detector, 
the data taking conditions, the kinematics of elastic J /tp production at HERA, and the 
event selection are described. The W dependence of the cross section cr^p^ji^p, the t 
distribution and the decay angular distributions are then presented. 



2 Experimental Conditions 
2.1 HERA 

During 1994 HERA operated with a proton beam energy of 820 GeV and a positron 
beam energy of 27.5 GeV. In the positron and proton beams 153 colliding bunches were 
stored together with 17 unpaired proton bunches and 15 unpaired positron bunches. The 
time between bunch crossings was 96 ns. The typical instantaneous luminosity was 1.5 x 

^q30 g-l_ 



2.2 The ZEUS Detector 

The main ZEUS detector components used in this analysis are outlined below. A detailed 
description of the ZEUS detector can be found elsewhere |jl5[. In the following the ZEUS 
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coordinate system will be used, the Z axis of which is coincident with the nominal proton 
beam axis, the X axis is horizontal and points towards the centre of HERA and the Y 
axis completes a right handed coordinate system. The origin of the coordinate system 
lies at the nominal interaction point. 

The momentum and trajectory of a charged particle were reconstructed using the Vertex 
Detector (VXD) ||] and the Central Tracking Detector (CTD) Q. The VXD and the 
CTD are cylindrical drift chambers which are placed in the solenoidal magnetic field of 
1.43 T produced by a thin superconducting solenoid. The CTD surrounds the VXD and 
covers the angular region 15° < 9 < 164° (where 9 is the polar angle with respect to the 
proton direction). 



The high resolution uranium-scintillator calorimeter CAL |T8| surrounding the coil is 
divided into three parts, the forward calorimeter (FCAL), the barrel calorimeter (BCAL) 
and the rear calorimeter (RCAL), which cover polar angles from 2.6° to 36.7°, 36.7° 
to 129.1°, and 129.1° to 176.2°, respectively. Each part consists of towers which are 
longitudinally subdivided into electromagnetic (EMC) and hadronic (HAC) readout cells. 

The proton remnant tagger (PRT), a set of scintillation counters surrounding the beam 
pipe at small forward angles, serves to tag events with proton dissociation. It is situated 
at Z = 500 cm and covers the angular range from 6 to 26 mrad. 



The muon detectors |[T^, situated outside the calorimeter, consist of limited streamer 
tubes (LST) placed both inside and outside the magnetised iron yoke. The inner chambers 
(BMUI and RMUI) were used to tag the muons from the J/ip. The BMUI and the RMUI 
cover the polar angles between 34° < 9 < 135° and 134° < 9 < 171°, respectively. 

Proton-gas events occuring upstream of the nominal interaction point are out of time 
with respect to the e~^p interactions and were rejected by timing measurements made by 
the scintillation counter arrays Veto Wall, C5 and SRTD situated along the beam line at 
Z = —730 cm, Z = —315 cm, and Z = —150 cm respectively. 

The luminosity was determined from the rate of the Bethe-Heitler process e~^p e'^'jp 
where the photon was measured by the LUMI calorimeter located in the HERA tunnel 
at Z = —107 m [^. The luminosity was determined with a precision of 1.5% for the 



measurements presented below. 

3 Kinematics 

Figure |l|a shows a schematic diagram for the reaction: 

e+{k)p{P) e+{k')J/tlj{V)p{P'), (1) 

where each symbol in parentheses denotes the four-momentum of the corresponding par- 
ticle. 

The kinematics of the inclusive scattering of unpolarised positrons and protons are de- 
scribed by the positron-proton centre of mass energy squared (s) and any two of the 
following variables 
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• Q"^ = —q^ = —{k — k'Y, the negative four-momentum squared of the exchanged 
photon; 

• y = {q ■ P)/{k ■ P), the fraction of the positron energy transferred to the hadronic 
final state in the rest frame of the initial state proton; 

• W"^ = {q + = + 2?/(/f; • P) + Mp fa i/s, the centre of mass energy squared of 
the photon-proton system, where Mp is the proton mass. 

For a complete description of the exclusive reaction e^p — > e^J/ipp {J/ip ^ where 
denotes a pair of electrons or muons) the following additional variables are required 

• t = {P — P'Y, the four- momentum transfer squared at the proton vertex; 

• the angle between the J/ip production plane and the positron scattering plane in 
the photon-proton frame, $; 

• the polar and azimuthal angles, 9h and (ph-, of the decay leptons in the J/ip rest 
frame. 

In the present analysis, $ is not measured because events were selected in which the 
scattered positron was not detected. In such untagged photoproduction events the 
value ranges from the kinematic minimum = M|?/^/(l — y) 10~^° GeV^, where 

Me is the electron mass, to the value at which the scattered positron starts to be observed 
in the uranium calorimeter Q'^^.x ~ ^ GeV^, with a median of approximately 5 x 
10~^ GeV^. Since the typical is small, the photon-proton centre of mass energy can 
be expressed as 

^ 2{Ej/^ - pzj/^)Ep = AEpE.y, (2) 

where Ep and Ej/^ are the laboratory energies of the incoming proton and the J/ijj and 
Pzj/ip is the longitudinal momentum of the J/ip. The four-momentum transfer squared, 
t, at the proton vertex for Q"^ = is given by 

t=iq-Vr^-plj/^, (3) 

where Ptj/iP is the momentum of the J/ip transverse to the beam axis. Non-zero values 
of cause t to differ from —p^jj^ by less than . A correction is applied to the p^j/^ 



distribution to correct for this effect as described in section |9l3| 111. 



4 Trigger 



ZEUS uses a three-stage trigger system [|T^. The electron and muon pair triggers are 
outlined below, followed by a summary of trigger requirements common to both channels. 
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Electron Channel 



The First Level Trigger (FLT) required 1, 2 or 3 track segments to be found in the CTD, 
with at least one segment pointing to the interaction region. The sum of all the energy 
deposited in the EMC section of the calorimeter was required to exceed 0.66 GeV. In 
addition, either the total energy in the calorimeter had to be greater than 2 GeV or the 
total energy in FCAL (ignoring the cells closest to the beam pipe) had to be greater than 
2.5 GeV. 

The Second Level Trigger (SLT) required the total energy in the HAG section of the 
calorimeter to be less than 1 GeV and the total energy in the EMC section to be greater 
than 1.5 GeV. The ratio of HAG to EMC energy in RCAL and BCAL separately had to 
be less than 0.1 or the HAG energy had to be less than 0.2 GeV. 

The Third Level Trigger (TLT) matched tracks measured in the CTD to electromagnetic 
energy deposits in the calorimeter. A cluster of contiguous cells, each with an energy of 
at least 0.3 GeV, was defined as electromagnetic if more than 90% of the total cluster 
energy was contained in EMC cells. An electron candidate was defined as a track with 
momentum transverse to the beam direction in excess of 0.4 GeV passing within 30 cm 
of the centre of an electromagnetic cluster. At least two electron candidates of opposite 
charge were required. At the distance of closest approach the separation between the two 
tracks was required to be less than 7 cm. An event was kept if the invariant mass of any 
pair exceeded 2 GeV. 

Muon Channel 

At the FLT, track segments had to be found in the inner barrel muon chambers (BMUI) 
accompanied by a reconstructed energy deposition of at least 0.464 GeV in a GAL trigger 
tower. Note that on average a muon produces a visible signal of 0.8 GeV in a trigger 
tower. Alternatively, hits had to be found in the RMUI chambers accompanied by a 
reconstructed energy deposit of at least 0.464 GeV in an RCAL trigger tower []15|. At 
least one and no more than five track segments had to be found in the CTD, with at least 
one pointing to the interaction region. 

No requirements were imposed at the SLT. 

At the TLT a muon candidate was formed when a track found in the CTD matched a 
cluster of energy in the calorimeter consistent with the passage of a minimum ionising 
particle (m.i.p.) and a track in the inner muon chambers. An event containing a muon 
candidate for which 6 > 147° was accepted if the momentum exceeded 1 GeV. The 
transverse momentum of a muon candidate for which 20° < 6 < 147° was required to 
exceed 1 GeV. 

Common Requirements 

An event was rejected at the FLT if the time of arrival of any signal observed in the Veto 
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Wall, the C5 counter or the SRTD was inconsistent with the time of the bunch crossing. 
In order to increase the purity of the sample the sum of energy in the inner ring of FCAL 
was required to be less than 1.25 GeV. 

At the SLT, the total energy in the calorimeter {Etoi = ^iEi) and the Z component of 
the momentum (Epz = SjE'jCos^j) was calculated. The sums run over all calorimeter 
cells i for which the energy, Ei, deposited in the cell is above threshold and the polar angle 
at which the cell is found is denoted by 6',. Beam-gas events were rejected by exploiting 
the excellent time resolution of the calorimeter. In order to remove inclusive beam-gas 
background in time with the bunch crossing, an event was rejected if the ratio T^pz/ Exot 
was greater than 0.96. 

Finally, at the TLT, Etoi and T^pz were calculated again using the CAL energies re- 
constructed at the TLT, and an event was accepted if Exot — ^Pz < 100 GeV and 
^Pz/Etoi < 0.94. 

5 Offline Event Selection 

To be accepted an event was required to have exactly two tracks of opposite charge with 
pseudorapidity, 77, in the range |?7| < 1.7. Denoting the polar angle of a track by 0, r] is 
defined such that 77 = — In (tan(6'/2)). The two tracks were required to fit to a common 
vertex consistent with an ep interaction. The tracks had to match to clusters of energy in 
the calorimeter and events were rejected if more than 1 GeV was deposited in calorimeter 
cells not associated with either of the two tracks. As shown in equation (^, W'^ was 
determined from the measured Ej/^ — pzj/^i} of the decay leptons. The requirement that 
the value of W lie in the range < W < 140 GeV restricted the sample to a region of 
high acceptance. Selection criteria specific to the electron and muon channel are described 
below. 

Electron Channel 

The electron sample comes from an integrated luminosity of 2.70 ± 0.04 pb~^. The algo- 
rithm used to define the electron pair sample at the TLT was reapplied offline with the 
final detector calibrations. The transverse momentum threshold of each of the two op- 
positely charged tracks was increased to 0.8 GeV. In order to reduce contamination from 
misidentified pions, the energy of at least one of the electromagnetic clusters matched to 
the tracks by the TLT algorithm applied offline was required to be larger than 1 GeV. 

Figure 0a shows the mass distribution of the electron pair sample. A clear peak at 
the J/ip mass is observed. The signal region, 2.85 < Me+e- < 3.25 GeV, contains 392 
events. The cross sections and angular distributions presented below are obtained by 
calculating acceptances and background contributions for this range. The solid line shows 
an unbinned likelihood fit in which a Gaussian resolution function has been convoluted 
with a radiative J/ip mass spectrum and a polynomial describing the background. The 
mass estimated by the fit is 3.094 ± 0.003 GeV, the rms width is 33 ± 4 MeV, and the 
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number of events attributable to J / ip production estimated by the fit over the mass range 
2 < Me+e- < 4 GeV is 460 ± 25. 



Muon Channel 



The muon sample comes from an integrated luminosity of 1.87 ±0.03 pb~ . The momen- 
tum of each track was required to exceed 1 GeV. At least one of the two tracks had to 
match a m.i.p. cluster in the calorimeter and a track segment in the barrel or rear muon 
chambers. To remove cosmic ray contamination the calorimeter signals were required to 
be in time with the beam crossing and the distance between the two tracks must be less 
than 2 cm at their distance of closest approach to the beamline. To further reduce the 
cosmic ray background the tracks were required not to be collinear. This was achieved 
by calculating the cosine of the angle, Q, between the two tracks at the interaction point. 
An event was rejected if cosQ < —0.99. 

The mass distribution for the events passing the muon pair selection is shown in figure 
Pd. a clear peak over a flat background is observed. The signal region, 2.95 < M^+^- < 
3.25 GeV, contains 289 events. The cross sections and angular distributions presented be- 
low are obtained by calculating acceptances and background contributions for this range. 
An unbinned likelihood fit to the sum of a Gaussian signal plus a flat background gives a 
value of 3.086 ±0.003 GeV for the mass, 38 ±3 MeV for the rms width and 266 ± 17 for the 
number of events attributable to J/ip production in the mass range 2 < M^+^- < 4 GeV. 



6 Monte Carlo Simulation and Acceptance Calcula- 
tion 

The reaction e"*" p — >■ e"*" J/ip p (figure |I]a) was modelled using the DIPSI Monte Carlo 



program This Monte Carlo is based on the model of Ryskin p2[ in which it is 

assumed that the exchanged photon fluctuates into a cc pair which then interacts with a 
gluon ladder emitted by the incident proton. The events are generated with a cross section 
proportional to and with an exponential t distribution proportional to exp{—b\t\). 
Good agreement between the generated and observed distributions is obtained for 6 = 1 
and 6 = 4 GeV~^. In order to determine the systematic error on the acceptance 6 was 
varied in the range < 5 < 2. The acceptance was found to be insensitive to the variation 
of b in the range 3 < 6 < 5 GeV~^. 

Events were generated in the W range 20 < W < 210 GeV and between Qf^i^ 

= 4 GeV^. The centre of mass decay of the J/ip was generated with a (1 ± a cos^ 9h) 
distribution with a = 1. Varying the value of a from 1 to 0.4, corresponding to about 
one standard deviation variation around the measurement presented in section |9.4|, the 



acceptance grows by less than 10%. A systematic error due to this uncertainty is included 
in the total systematic error as described in section ^. The effects of positron initial and 
final state radiation and that of vacuum polarisation loops were neglected; the effects on 
the integrated cross section have been estimated to be smaller than 4% [|l|. 
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The events were then passed through a detailed simulation of the ZEUS detector and 
trigger. Parameterisations of noise distributions obtained from data taken with a random 
trigger were used to simulate the calorimeter noise contribution to the energy measure- 
ments. The simulated events were subjected to the same reconstruction and analysis pro- 
grams as the data. The distributions of the reconstructed kinematic quantities obtained 
using DIPSI are in good agreement with those from the data. The overall acceptance 
was obtained as the ratio of the number of accepted Monte Carlo events to the number 
generated in the selected kinematic range. The acceptance, calculated in this manner, 
accounts for the geometric acceptance, for the detector, trigger and reconstruction effi- 
ciencies, and for the detector resolution. Table |^ shows the acceptances in various W 
ranges determined for each decay mode. 



7 Background 

In addition to elastic J/ ip photoproduction, the following processes may contribute to the 
final sample: 

• The Bethe-Heitler process in which a lepton pair is produced by the fusion of a 
photon radiated by the positron with a photon radiated by the proton. This process 
was simulated using the LPAIR Monte Carlo ||2^ which was used to generate events 
in which the proton remains intact ('elastic' events) and events in which the proton 
dissociates ('dissociative' events). The size of the Bethe-Heitler contribution to the 
non-resonant background is shown in figure ^ where the t^l~ mass distributions are 



plotted. The QED cross section pj| for the elastic and dissociative Bethe-Heitler 
processes have been used to determine the normalisation of the appropriate LPAIR 
Monte Carlo sample. Figure ^ shows that the Bethe-Heitler process saturates the 
non-resonant background in the muon channel and is the dominant source of non- 
resonant background in the electron channel. The calculated background due to 
the Bethe-Heitler process in the signal region is 38 ± 1 for the electron channel and 
23 ± 1 for the muon channel. 

Pions misidentified as electrons in the electron sample. For e^e~ masses larger than 
2.5 GeV the Bethe-Heitler contribution saturates the non-resonant background. The 
residual contribution of misidentified pions in the final sample was shown to be less 
than 1.5% by studying the distribution of dE/dX obtained using the pulse height 
information from the CTD. No subtraction has been made for pion misidentification. 
A systematic error of —1.5% attributed to the uncertainty in the pion contamination 
was included in the final systematic error. 

J ftp produced via the production and decay of ip'. The only i/j' decay mode giving 
a significant contribution to the J/tp signal is tp' J/ipir^TC^. 

Proton dissociative J/ip production (figure |T|b). The EPSOFT Monte Carlo was 
used to simulate this process. EPSOFT is based on the assumption that the diffrac- 
tive cross section is of the form da/d\t\dMff oc e"'''^'*'/^^ where Mn is the mass 
of the dissociative system. The simulation of the hadronisation of the dissociative 
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system includes a parameterisation of the resonance spectrum. To cross-check the 
results the generator PYTHIA was also used which contains a different param- 
eterisation of the resonance spectrum. 



After the subtraction of the Bethe-Heitler contribution, the production of J/ip mesons 
via the decay of the i})' and proton dissociative J/ip production are the only significant 
sources of background and will be discussed separately below. 

The ip' contribution was determined using a sample of events in which the ip' decayed 
to a muon pair (branching ratio Bi = (0.77 ± 0.17)% |2^). This sample was obtained 
using the same cuts as those used to isolate the J/ip ^ jj,^ ji^ sample (see section |]). A 
signal of A^^i = 7 ± 4 events was found at the ip' mass in a sample for which the integrated 
luminosity, £i, was 2.70 ± 0.04 pb~^. The corresponding acceptance, Ai, computed with 
DIPSI, was Ai = 0.35. The number of events from ip' production entering the elastic 
J/ip ^ fi^fi^ sample via the decay ip' — ^ J/ipir^n^ was estimated using the formula 

Nc = j-^^^c^cBcB, (4) 

where B = (6.01 ± 0.19)% is the branching ratio for the decay J/ip /W^/U^, Be is that 
for the decay ip ipn^n^ {Be = (18.4 ±2.7)%) Cc is the luminosity from which the 
muon sample defined in section ^ was drawn [Cc = 1-9 pb~^) and A(j is the acceptance 
for the process e'^p — > e^ip'p {ip' — » fi~^ fi~7T^TT^), using DIPSI Aq = 0.28. The formula (^ 
leads to a ip' contamination of (2.3 ± 1.4)%. This result was cross-checked by selecting 
events in which the ip' decayed into /i"'"/i~7r+7r~. In this case 7 ± 3 events were found at 
the Ip' mass and a contamination of (3.4 ± 1.4)% was estimated. The two results may be 
combined to give a final estimate of the ip' contamination of (3± 1)%. This contamination 
was subtracted from both the electron and muon sample. 

The proton dissociative process is characterised by a cross section of the form 
da e"^''!*! 

d\t\dMi; ^ 1^- 

In order to estimate the value of 6^, dissociative events were selected in which the J /ip was 
accompanied by an energy deposit in the inner ring of FCAL or in the PRT. The value 
hd = 1 GeV~^ was found to give the best description of the Ptj/^ distribution of the 
PRT tagged sample. The systematic error in the dissociative contribution caused by the 
uncertainty in 6^ was estimated by varying bd in the range 0.4 < bd < 2 GeV~^. This 
assumption is consistent with the result bd = 1.6 ± 0.3 ±0.1 GeV~^ reported by the HI 
collaboration [T^. The value /? = 2.25 was used as the central value in the simulation of 



the Mat distribution and (3 varied in the range 2 < /3 < 2.5 to estimate the systematic 
error. This assumption is consistent with the result /3 = 2.20 ± 0.03 recently obtained at 
Fermilab for the diffractive dissociation of the proton in pp collisions The mass of 
the nucleonic system was generated in the range (1.25 GeV^) < < 0.1 W"^. 

The proton dissociative contribution to the electron sample was determined by selecting 
a sample, D^, for which the requirement that Exot — Ej/^ < 1 GeV was replaced by the 
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three cuts Ep > 1 GeV, Eb < 1 GeV and Er < 1 GeV. Ep, Eb and Er were calculated 
by summing the energy in the FCAL, BCAL and RCAL respectively. The calorimeter 
cells associated with the electron candidates were excluded from these sums. The cut on 
Ep selects dissociative events in which energy is deposited in the proton direction, while 
the cut on Er ensures that events in which the scattered positron is detected in RCAL 
do not enter the sample. The cut on Eb ensures that inelastic events depositing energy 
in BCAL also do not enter the sample. The proton dissociative sample, Dg, was further 
examined by studying the distribution of the energy weighted pseudorapidity defined by 



TiiEiTji , . 

= ^ p , (6) 



where Ei is the energy of a calorimeter cell and r]i is the pseudorapidity of the cell and the 
sum nms over all cells containing more than 200 MeV but excluding those matched to the 
tracks forming the J/ip candidate. The distribution of fjc for dissociative events, simulated 
using the EPSOFT Monte Carlo, is strongly peaked at r^c > 2. In the sample D^. there are 
2 events for which fjc > 2. The ratio of the number of EPSOFT events passing the elastic 
cuts to the number with Ep > 1 GcV, Eb < 1 GcV, Er < 1 GcV and fjc > 2 was 58. This 
leads to a dissociative contribution to the clastic J/ip to electron sample of (331^2^6^18)%- 
The first error is statistical and the second error is the systematic error resulting from 
the allowed variation of P in the Monte Carlo generation of dissociative events. When 
the calculation is repeated with EPSOFT replaced by PYTHIA the result differs by - 
18% from that reported above. The third error quoted in the dissociative contribution 
reflects this uncertainty in the simulation of the dissociative flnal state. The change in the 
dissociative contribution obtained when bd was varied in the range 0.4 < bd < 2 GeV~^ 
was found to be negligible. 

The same procedure was applied to the muon sample with the only difference that the cut 
on f]c was not applied. The proton dissociative sample obtained contained 7 events and 
the ratio of the number of EPSOFT events passing the elastic cuts to the number with 
Ep > 1 GeV, Eb < 1 GeV, Er < 1 GeV was 11. This leads to a dissociative contribution 
of (29±llt^ 

Independent estimates of the dissociative contribution were made using dissociative events 

tagged by the PRT. EPSOFT was used to estimate the fraction of untagged dissociative 
events in the elastic sample since it was found that PYTHIA gives a poor description of the 
multiplicity distribution observed in the PRT. The dissociative contamination estimated 
in this way was (34 ± 8)% for the electron channel and (27 ± 8)% for the muon channel. 
The errors quoted are statistical only. 

The four independent results were combined to give a final estimate of the dissociative 
contribution of (30 ± 5tl tio)%- 



8 Systematic Errors 

Several factors contribute to the systematic errors in the elastic J/ip cross section measure- 
ment. In the following they are divided in two categories: decay channel specific errors 
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and common systematic errors. The first category contains systematic errors specific to 
the electron or muon decay channel, while the second contains systematic errors common 
to both decay channels. Table ^ summarises all these systematic errors. 



Decay channel specific errors: 



Trigger. For the electron channel, the dominant systematic error due to the FLT 
acceptance is given by the requirement Exot > 2 GeV. At the SLT the dominant 
systematic error is contributed by the simulation of the calorimeter noise. For the 
muon channel, the dominant systematic error is contributed by the uncertainties in 
the simulation of the trigger threshold and the CTD-FLT track reconstruction. No 
systematic error in either channel is attributed to the TLT acceptance since all cuts 
are superseded by more stringent requirements offline. 

Event selection: In this class we include the systematic errors due to uncertainties 
in the measurement of momentum, transverse momentum, \r]\ and the choice of the 
mass window. For the electron channel uncertainties in the cuts used to define an 
electron cluster also contribute. For the muon channel this class also contains the 
uncertainties coming from the collinearity cut. Each cut was varied within a range 
determined by the resolution of the quantity in question and the changes induced 
in the results were taken as an estimate of the corresponding systematic error. The 
different systematic errors were summed in quadrature. 

Pion misidentification: This class applies to the electron channel only; the method 
used to determine the systematic error was described in section |^. 

Muon chamber efficiency: The systematic error attributed to errors in the muon 
chamber reconstruction efficiency was estimated by using cosmic ray events. 



Branching ratio: The error on the branching ratio J/ip i^i as quoted in [E6|. 



Common systematic errors: 



• Acceptance: The uncertainty in the acceptance was estimated by varying the pa- 
rameters b and S as described in section ^. 

• Elastic definition: The systematic uncertainty contributed by the criterion used to 
classify an event as elastic was estimated by changing the elastic definition: Exot — 
Ej/^ < 1 GeV to Erot - Ej/^ < 0.7 GeV and to Etoi - Ej/^ < 1.3 GeV. 

• Radiative corrections: The effects of positron initial and final state radiation and 
that of vacuum polarisation loops were neglected; the effects on the integrated cross 
section have been estimated to be smaller than 4% [|I|. We take 4% as an estimate 
of the systematic error attributable to this source. 

• Helicity distribution: The centre of mass decay of the J/ip was generated with a 
(1 + acos^Oh) distribution. The systematic error was evaluated by varying the 
value of a from 1 to 0.4. 
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Mat distribution in proton dissociation: As explained in section ^ this is obtained 
by changing the parameter f3 in the range 2 < (3 < 2.5. 

Model of dissociation: The dependence on the modelhng of the hadronic final state 
in proton dissociation was obtained by comparing the contamination obtained using 
PYTHIA with that obtained using EPSOFT (see section 0). 

ijj' contamination: As explained in section |^ the systematic error on the ijj' contri- 
bution is 1%. 

Luminosity: As indicated in section P]2| the uncertainty of the luminosity determi- 
nation is 1.5%. 



9 Results 



9.1 Integrated Cross Sections 

The cross section for elastic J/ip electroproduction is given by 

_ NEvt ,„s 

where C is the integrated luminosity, A is the acceptance, B is the branching ratio for 
J/ip to decay into electron or muon pairs |^ and N^vt is the number of signal events 
after background subtraction. NEvt and A were determined in the signal regions defined 
for the electron and muon channels in section ^. In the range 40 < W < 140 GeV and 
for Qfnin < < 4: GeV^ the J/ip electroproduction cross section is 

aep^ej/ijp = 5.37 ± 0.30(stat.)^°;gg(syst.);^o^'^(model) nb, (8) 

using the electron sample and 

(Tep-,ej/^pp = 5.04 ± 0.32(stat.)^Q;7g(syst.);!;Q^°(model) nb, (9) 

using the muon sample. The model error quoted above is due to the difference between 
the value of the dissociative contribution estimated using EPSOFT and using PYTHIA. 
In the systematic error we have summed in quadrature all the decay-channel-specific 
errors and the common systematic errors. The electron and muon cross section results 
are compatible with each other and with previous measurements in the same W range 

mmm. 
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9.2 Photoproduction Cross Section 



The photoproduction cross section is related to the ep cross section by p8 



my,Q')dydQ^ " $t ' ^ ' 

where cr^p_j/^p is the mean cross section in a range of W and $t is the effective flux of 
virtual photons accompanying the positron. The integrals run over the full range of 
and from i/min = W^^^/s to i/max = W^^^/s where Wmin and Wmax are the minimum 
and maximum values of W respectively. The photoproduction cross section has been 
determined in four W bins. The results for each of the lepton decay modes and the 
combined results are reported in table [1|. The procedure described in section was used 



to calculate the errors on the cross sections presented in table |I]. For the combined results 
the following procedure was used. The weighted mean cross section was calculated; the 
weights being obtained by summing the statistical and decay channel specific errors in 
quadrature. The first error reported on the combined results in table [l| is the error on the 
weighted mean, the second error is the sum of the common systematic errors added in 
quadrature. The third error reported on the combined results in table |l] is the systematic 
error associated with the model of diffraction. The combined results are shown in figure |^ 
where a^p^j/^p is plotted as a function of W. The points are plotted at the mean values 
of W reported in table |1|. A clear growth of a^p^j/^p with W is observed over the W 
range covered by this experiment. 

The ZEUS data in the range 40 < W < 140 GeV were fit to the form a^p^j/^p oc 
with the result S = 0.92±0.14 (stat.)±0.10 (syst.). The systematic error was obtained as 
follows. For each source of systematic error in turn the cross sections were displaced from 
their central values, the fit was performed and the value Ssi recorded. The systematic 



error on 6 was taken to be yJ2iiS — Sgi) ■ The result of the fit is shown in figure ||a. 
This value of 6 disfavours that expected in the Donnachie-Landshoff model |^ (the soft 
pomeron model) in which 6 is expected to take the value 6 = 0.22 in this W range. The 
curve corresponding to the soft pomeron model is shown in figure 3a as a dotted line 
arbitrarily normalised to the second ZEUS data point. 

It is interesting to compare the ratio, i?(^),ofthe cross section for elastic J/ ip photopro- 
duction to the cross section for elastic p production as a function of W. At ~ 12 GeV 
R{^) = (1.21 ±0.20) X 10-3 while at ~ 15 GeV i?(^) = (1.67 ±0.23) x 10"^ 
[p!0| , [Tl|, ^ . The results presented in the present paper may be combined with those pre- 
sented in reference to determine that i?(^) = (2.94 ± 0.74) x 10"^ at ~ 70 GeV 

showing that R[^-y-) rises with W . These values are to be compared with R['-^) = | 
expected on the basis of the quark charges and a flavour independent production mecha- 
nism. 

The data are replotted in figure §3 together with other measurements of elastic J/ip pho- 
toproduction. The results of two pomeron models |32| are shown in figure In 



the model of reference |^ the effective pomeron intercept is assumed to depend upon 
Q^HKK = cM^ + Q^, where Mc is the mass of the charm quark and the constant c ~ 1. 
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The model of reference assumes a fixed pomeron intercept but includes both a scale 
dependent pomeron couphng and a mass threshold function. Both models give a good 
description of the data. 

Attempts have been made to describe elastic J /ip production in perturbative QCD, pQCD. 
In the approach of Ryskin |^ the pomeron is described as a gluon ladder evaluated in 



the leading logarithm approximation. In this model the cross section is proportional to 
[asXg{x,(f)f' , where is the strong coupling constant (assumed fixed and set equal to 
0.25) and xg{x, (f) is the gluon momentum density in the proton. The quantities x and 
(f are given by 

^= m ^= 4 ^^^^ 

and give the effective momentum fraction and scale at which the gluon density is probed 
respectively. In the present case both and \t\ are negligible in comparison to Mj^^. 
For elastic J/ip photoproduction takes a value of approximately 2.5 GeV^ while 
the measurements presented here are sensitive to values of x in the range 0.4 x 10^^ < 
X < 6 X 10~^ 0]. If a gluon distribution of the form xg (x, Q^) oc x~^ is assumed then 
the W dependence of a^p^j/^p may be written o^p^jj^p oc VT^'^. The value of b reported 
above gives A = 0.23 ± 0.04 ± 0.03. This is consistent with our measurement of the gluon 
distributions based on an analysis of the scaling violations of F2 extrapolated back to 
Q2 = 2.5 GeV^ p3[. 



Figure shows the results of the pQCD calculation of a^p^j/^p presented in |3^ which 



extends the Ryskin model beyond leading order and includes the effects of the relativistic 
motion of the c and c within the J / ip and the rescattering of the cc pair on the proton. 



Good agreement with the data is obtained using the MRS-A' ||3^ parton distributions. 
Other choices of parton distributions compatible with HERA measurements of F2 also give 
an acceptable description of the W dependence of a^p^j/^pp over the range 40 < < 140 
GeV. 



9.3 Differential Cross Sections 

Figure Ija shows the differential photoproduction cross section da/dp^jj^ for the full W 
range (40 <W < 140 GeV). The results from the electron and muon samples have been 



combined using the procedure described in section |9.2| . The contribution from proton 
dissociative J/ijj production and the Bethe-Heitler process have been subtracted bin by 
bin. The cross section exhibits the exponential fall characteristic of diffractive processes. 
A binned likelihood fit to the form 

= Ae-'^^"'^^/^ (12) 

was performed in which the function in equation [l^ was integrated and compared with 
the measured cross section bin by bin. Fitting over the range p^j/^ < 1 GeV^ gives the 
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result 



bp^ = 4.3 ± OAt^oi GeV-2. (13) 

The differential cross section da/d\t\ may be obtained by dividing da/dp'^j^^ bin by bin 
by a factor which corrects for the small of the photon. Figure ^ shows the correction 
factor, F, which is slowly varying and close to 1 for \t\ < 1 GeV^. The differential cross 
section da/d\t\ obtained in this way is plotted in figure Again, the cross section 
exhibits an exponential fall and a binned likelihood fit to the form 

4^ = ^e-^l*l (14) 
d\t\ ^ ' 



was performed in which the function in equation [Tj was integrated and compared with 
the measured cross section bin by bin. Fitting over the range |t| < 1 GeV^ gives the result 



6 = 4.6 ± 0.4+°;^ GeV-^. 



(15) 



The systematic error contains the contribution coming from the uncertainty in the cor- 
rection factor F. The fit for b was repeated for \t\ < 0.8 GeV^ and |t| < 1.2 GeV^. 
The small changes in b obtained are included in the systematic error quoted in equation 
|T^. The size of the statistical and systematic errors on the parameter b prevents us from 
investigating the dependence of 6 on using the data presented here. The slope is in 
agreement with the result obtained by the HI collaboration |1^, |13| in the same W range. 
We have previously determined the parameter b in elastic p, u and photoproduction 
to be 9.8 ± 0.8 stat. ± 1.1 syst. GeV-^ |l|, 10.0 ± 1.2 stat. ± 1.3 syst. GeV-^ |] and 
7.3 ± 1.0 stat. ± 0.8 syst. GeV^^ respectively. In geometrical models of vector meson 
production these results may be interpreted as indicating that the radius of the J/ip is 
smaller than that of the p, uj and 0. When the parameter b is measured in exclusive p 
production in deep inelastic scattering for Q"^ values in the range 7<,Q^ <,25 GeV^ a value 
of 5.l;*^o]9 ^ 1 GeV~^ is obtained which is significantly smaller than the slope obtained in 
elastic p photoproduction . Thus, in exclusive p production b falls as is raised from 
reaching a value of 5.11^'^— 0.9 ± 1.0 GeV~^ at of order 10 GeV^ comparable to that 
reported here for J /if) photoproduction where the hard scale in the scattering process may 



be set by Mj 



9.4 Decay Angular Distributions 



The J/ ip decay angular distributions can be used to determine elements of the J/ if) spin- 
density matrix . In the s-channel helicity frame the J / is at rest and the quantisation 
axis is taken to lie along the J/ip direction in the photon-proton centre of mass system. 
The decay angular distribution is a function of 6^ and (ph, the polar and azimuthal angles 
of the positive lepton in the helicity frame. The angular distributions can be shown to be 
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dN 



N d cos 9h 



1 + r°^ 

' 00 



3C^) cos^^, 



(16) 
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^ ^ [l + r?i,cos20j. (17) 



N d(ph 271 

In the present experiment ^ GeV^ so that the J/tp is expected to be produced 
predominantly by transverse photons. If s-channel helicity is conserved (SCHC) then 
the spin density matrix parameters Tqq and r^i^^ should be close to zero. Under the 
assumption of SCHC, Tqq can be related to the ratio of the photoproduction cross sections 
for longitudinal and transverse photons 

^-^T^ ■ (18) 

where e is the virtual photon polarisation, i.e. the ratio of the flux of longitudinally 
polarised photons to the flux of transversely polarised photons. The mean value of e over 
the kinematic range of and y sampled by the present experiment is e = 1.043. 

The decay angular distributions are presented in figure ^. No subtraction of the disso- 
ciative contribution has been made for the distributions presented in figure ^ since it is 
assumed that the elastic and dissociative processes have the same angular dependence. 
The distribution of the polar angle of the positive lepton is shown in figure |^. A fit has 
been used to determine Tqq. The result Tqq = —0.01 ±0.09 (which gives R = —0.01 ±0.09) 
is consistent with SCHC. The distribution of the azimuthal angle of the positive lepton 
is shown in figure §b. The distribution is fiat and a fit has been made to determine r^^i. 
Again, the result r°i^ = -0.08 ± 0.07 is consistent with SCHC. 



10 Summary 

The cross section for elastic J ftp photoproduction has been measured using the ZEUS 
detector at HERA. A significant rise in the cross section with W has been observed for 
W in the range AO < W < 140 GeV. The rise in the cross section with W may be 
parameterised by cx-yp^j/^p oc with 5 = 0.92 ± 0.14(stat.) ± 0.10(syst.). The measured 
value of S is inconsistent with the soft pomeron model. Models based on the vector 
dominance model plus the exchange of a pomeron can be made to describe the data if 
the effective pomeron intercept, or the effective pomeron coupling is assumed to depend 
on the hard scale in the process. QCD based models, which describe the process in terms 
of the exchange of a gluon ladder evaluated at leading order or beyond leading order, are 
consistent with the data. 

The differential cross section da/d\t\ has been measured and falls exponentially with \t\. 
The slope of the exponential has been measured to be 4.6 ± 0.4lQ;g GeV~^ in the range 
|t| < 1 GeV^. In geometrical models of vector meson production these results may be 
interpreted as indicating that the radius of the J/ip is smaller than that of the p, uj and 
(j) as measured in photoproduction. 

The decay angular distributions are consistent with s-channel helicity conservation. 
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W Range 
(GeV) 


Mode 


^ Dig 


A 


(Jcr, J Iti-i'n ( lib ) 




(J^r, i. T /th'Tt I lib ) 


(^-fp^j/ipp (lib) 
combined 


40-60 
{W) = 49.8 ±0.8 




84 ± 10 


0.28 


1 23 ± 14+°-^I +°n'' 

J..^KJ _1_ »-'.-^t:_o.21 — 


0.0411 


29 9 ± 3 4+ + r 


30 4 ± 3 4+^-^ + 




48 ± 7 


0.23 


1.28 ±0.19+^-^^ 


0.0411 


31.1±4.6+H + F 


60-80 
(W) = 71.2 ±0.7 




98 ± 11 


0.33 


1.24±0.13+^i^ 


0.0266 


46.6 ± 4.9+ + 


42.9 ± 4.5+H + f;^ 




61 ± 8 


0.35 


1.05±0.14+°13 +0^11 


0.0266 


39.5 ± 5.3+^^ + 


80-100 
{W) = 89.6 ±0.7 




92 ± 10 


0.32 


1.19±0.13+^-}I 


0.0189 


63.0 ± 6.9V.'' +0 ^ 


57.7± 5.8+«Q +^n^ 




70 ± 9 


0.42 


1.01 ± 0.12+° }1 


0.0189 


53.4 ± 6.3l?:l 1 


100-140 
{W) = 121 ± 1 




81 ± 9 


0.21 


1.59 ± 0.18+!]-i +V' 


0.0251 


63.3 ± 7.2+?nV V 


66.5 ± 6.8+^i +^n^ 




87± 10 


0.30 


1.74±0.20l°:i IV^ 


0.0251 


69.3±8.0l?iV-1) 



Table 1: The results for the integrated J/'^ photoproduction cross section as a function of W . Nsig is the number of events after 
subtraction of the Bethe-Heitler contribution and A is the acceptance. The photon flux $t is calculated as described in the text 
and used to calculate the 7p cross section, ci^p^^j/^p, from the ep cross section, Cep^ej/^p- Cross sections for the individual channels 
are quoted with the first error being statistical and the second systematic. The third error is the error attributed to the model of 
proton dissociation used for background subtraction and is described in the text. The combined electron and muon results have been 
obtained by averaging as described in the text. Here the first error contains the combined statistical and decay channel specific errors 
while the second contains all sources of common systematic error. The error attributed to the model of proton dissociation is the 
third error. 



Breeikdown of Contributions to the Systematic Error 

Values are quoted in percent 
Decay Channel Specific Systematic Errors 




Electron Channel 


Muon Channel 


W bin (GeV) 


40-60 


60-80 


80-100 


100-140 


40-60 


60-80 


80-100 


100-140 


Trigger 


+7 

-7 


+7 

-7 


+7 

-7 


+7 
-7 


+5 
-5 


+5 
-5 


+5 
-5 


+5 
-5 


Event selection 


+5.7 
-4.6 


+4.6 
-4.9 


+3.0 
-4.0 


+8.8 
-3.8 


+■2.0 
-6.3 


+5.6 
-4.1 


+7.1 
-0.0 


+7.0 
-2.0 


Pion misidentification 


+ 
-1.5 


+ 
-1.5 


+ 
-1.5 


+ 
-1.5 










Muon chamber 
efficiency 










+2 
-2 


+2 
-2 


+2 
-2 


+2 
-2 


Branching ratio 


+:i.-2 

-3.2 


+:i.-2 

-3.2 


+3.2 
-3.2 


+3.2 
-3.2 


+3.2 
-3.2 


+3.2 
-3.2 


+-i.2 
-3.2 


-3.2 


Subtotal 


+9.6 
-9.1 


+9.0 
-9.2 


+8.3 
-8.8 


+11.7 
-8.7 


+6.6 
-8.9 


+8.4 
-7.5 


+9.5 
-6.3 


+9.4 
-6.6 


Common Systematic Errors 


W bin (GeV) 


40-60 


60-80 


80-100 


100-140 


Acceptance 


+3 
-3 


Elastic definition 


+1 
-3 


Radiative corrections 


+4 
-4 


Helicity distribution 


+ 
-10 


+0 

-8 


+0 
-6 


+ 
-10 


Proton dissociation 


+6 
-7 


Model of dissociation 


+10 
-0 


tp contamination 


+1 
-1 


Luminosity 


+1.5 
-1.5 


Total Systematic Errors 




Electron Channel 


Muon Channel 


W bin (GeV) 


40-60 


60-80 


80-100 


100-140 


40-60 


60-80 


80-100 


100-140 


Total 


+16.0 
-16.4 


+15.7 
-15.3 


+15.3 
-14.1 


+17.4 
-16.2 


+14.5 
-16.3 


+15.4 
-14.4 


+16.0 
-12.7 


+15.9 
-15.2 



Table 2: The contributions to the systematic errors on the J/ip photoproduction cross 
section. The contributions to the systematic error are divided into Decay Channel Specific 
Systematic Errors and Common Systematic Errors as described in section 8. 
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Figure 1: Schematic diagrams for diffractive J/tl) electroproduction. (a) The mechanism 

for elastic vector meson production, (b) Proton dissociative J / i/j photoproduction where 
the proton dissociates into a hadronic system of invariant mass Mjv- 
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Figure 2: (a) The mass distribution of the events in the electron pair sample. A clear 
peak at the J/ if} mass is observed. The solid line shows the result of a fit in which a 
Gaussian resolution function has been convoluted with a radiative J/%1) mass spectrum 
and added to a polynomial background, (b) The mass distribution for events in the muon 
pair sample. The solid line shows the result of a fit in which a Gaussian resolution function 
has been added to a flat background function. For both the electron and muon channels 
the contribution of events from the Bethe-Heitler process is shown as the hatched area. 
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Figure 3: The elastic J/ip photoproduction cross section as a function of W. (a) Shows 
the results of this analysis. The error bars represent the statistical and systematic errors 
added in quadrature. The solid line shows the result of the fit to the data using the 
expression a^p^j/^p oc . As described in the text the value 5 = 0.92 ± 0.14 ± 0.10 
was obtained. The dashed line shows the prediction of a soft pomeron model in 
which 6 ~ 0.22. (b) The results of this analysis (solid circles) are compared to data from 
HI, ZEUS and the results of lower energy measurements |]^, |TU]. The result of a pQCD 
calculation |3^ in which the MRS-A' parton distributions have been used is shown 
as the sohd line. The result of the calculation presented in ||3l| is shown by the dotted 



line. The result of the calculation presented in p2| is shown by the long dash dotted line. 
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Figure 4: The distribution of transverse momentum squared for J/ip produced in the 
reaction 7p — > J/t/jp in the kinematic range AO <W < 140 GeV. (a) The differential cross 
section da/dp^j^^. The data are shown as the points and the result of the exponential 
fit in the range < 1 GeV^ is shown as the solid line, (b) The correction factor, F, 
required to obtain the \t\ distribution from the p%ji^ distribution by accounting for the 
of the photon, (c) The differential cross section da/d\t\. The result of the exponential 
fit in the range |t| < 1 GeV^ is shown as the sohd fine. In (a) and (c) the inner error 
bars represent the statistical and decay-channel-specific errors added in quadrature, the 
outer ones statistical, decay-channel-specific errors and common systematic errors added 
in quadrature. 
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Figure 5: Acceptance corrected decay angular distributions for the J/ip in the reaction 
ep — > eJ/ipp in the kinematic range 40 < < 140 GeV. No subtraction of the proton 
dissociative contribution to the sample has been made for the data presented in this 
figure since the angular dependence of the proton dissociative and elastic J/ production 
is assumed to be the same. The curves arc the results of the fits described in the text. The 
error bars represent the statistical, decay-channel-specific errors and common systematic 
errors added in quadrature. 
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